Properties of the Two-Dimensional Electron Gas Close to the Fermi-Liquid Quantum 
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The rearrangement of single-particle degrees of freedom of a dilute two-dimensional electron gas 
in the vicinity of the quantum critical point is examined within a microscopic approach. It is shown 
that just beyond the critical point, the Landau state undergoes self-consistent rearrangement of 
the quasiparticle spectrum and momentum distribution. At very low temperatures, there emerges a 
multi-connected quasiparticle momentum distribution. With increasing temperature, two crossovers 
occur: a fermion condensate appears in the first and disappears in the second, giving way to universal 
non-Fermi-liquid behavior. Manifestations of these crossovers in thermodynamic properties of the 
electron gas are studied and characterized. The four quasiparticle phases predicted to exist in the 
vicinity of the critical point are collected in a schematic phase diagram. 

PACS numbers: 71.10.Hf, 71.27.-|-a, 71.10.Ay 



INTRODUCTION 



The two-dimensional (2D) electron gas is a widely used 
physical model for the electron system of the inversion 
layer in high-quality Si-MOS devices. The density of 
that system can be varied over several order of magni- 
tude upon varying the external electric field orthogonal 
to the layer. The potential energy of the 2D electron gas 
is proportional to the inverse of the average distance be- 
tween electrons, while its kinetic energy is proportional 
to the inverse square of this distance. Thus, with decreas- 
ing density, the electron gas passes from a regime of weak 
correlations to one of strong correlations. Although ex- 
perimental study of the 2D electron gas began almost 
40 years agoji it is only recently that experiments^ 
have reached a degree of accuracy such that the explo- 
sive growth of the effective mass M* of Landau quasi- 
particles is observed as the density approaches a crit- 
ical value poo — 8 X lO^^cm"^. Evidence for diver- 
gence of the quasiparticle effective mass is also found 
in 2D liquid '^He (see Refs. UEH) and in some heavy- 
fermion materials i^i^ d^d^d^ At zero temperature, the 
point X = Xoo of an external parameter x (density, pres- 
sure, magnetic field) at which the effective mass is pre- 
dicted or found to diverge is known as the quantum crit- 
ical point (QCP)iiiiii^ Properties of strongly correlated 
Fermi system in the vicinity of the QCP exhibit non- 
Fermi liquid (NFL) behavior, i.e. they cannot be de- 
scribed within the Landau FL theory, as ordinarily in- 
terpreted. It is often supposed (see e.g. Refs. [l3fl5lfT6l ) 
that the collective degrees of freedom are the dramatis 
personae near the QCP; indeed, that they are respon- 
sible for the occurrence of the QCP itself and for the 
properties of the system in its vicinity. We follow a dif- 
ferent strategy, exploring instead the possibility that the 



single-particle (sp) degrees of freedom are the main char- 
acters of the drama. In this paper, we will show that the 
QCP is inherent in Landau FL theory, and that a tool 
exists within this theory for describing the system in the 
immediate vicinity of the QCP and on both sides of it. 

We focus on the 2D electron gas, for which a micro- 
scopic, ab initio non-pert urbative approach was devel- 
oped in Refs. IT8lfT9ll20l . Applying this method, it was 
found that at T = the effective mass M* of the 
2D electron gas diverges as the dimensionless parame- 
ter rs = ^/2Me^ /pf (where is the Fermi momentum) 
approaches a critical value ~ 7. This result for the 
critical spacing parameter is in reasonable agreement 
with the experimental value ~ 8 — 9 cited in Refs. ^H^. 

Fig. 1 presents results^ of evaluation of the single- 
particle energy e(p) (measured from the chemical poten- 
tial fi) , for the 2D electron gas at a set of rg values near 
the critical point r^. As seen in this figure, the spec- 
trum acquires an inflection point at the critical dilution, 
i.e., e(p, r^) (x (p — pfY, a situation first analyzed in 
Ref. [2l|. Beyond this point, meaning > r^, the con- 
ventional Landau state, characterized at zero tempera- 
ture by a Fermi-step quasiparticle momentum distribu- 
tion nFL(p) = 9{pf —p), becomes unstable against spon- 
taneous creation of quasiparticle-quasihole pairs. We call 
attention here to the fact that at T = a necessary condi- 
tion for stability of the conventional Landau state is pos- 
itivity of the variation of the ground-state energy with 
respect to any admissible variation of the quasiparticle 
momentum distribution from n-pi^{p). More explicitly, it 
is necessary that 



SE = J e(p, [n]) 5n{p)dv > (1) 
holds for n{p) = npi^ijj) and variations satisfying the 



2 



0.4 



0.2 



0.0 



-0.2 



i i i i 1 i 
2D electron gas 


1 r-» — ji T 

,''//. 




,' / / / 


r=8 


// / / ' 




■ ',' / ^ 


: .'\___ 


' >/ / 






- . - ' ' ' '^=5 





0.5 



1.0 
P/Pp 



1.5 



FIG. 1: Single-particle spectrum e(p) of the homogeneous 
two-dimensional electron gas in units of e% — p%./2M, evalu- 
ated at T = for different values of r . . 



particle-number constraint 

5n{p) dv — , 



(2) 



where dv — 2d'^p/ {2tt)^ . This stability condition ([T]) is 
met by ripL provided the corresponding spectrum e(p) has 
a single zero at p = pp. But as seen from Fig. 1, just be- 
yond the critical point the spectrum takes on a cubic- 
like behavior with three zeros near the Fermi surface. (It 
is worth noting that a similar shape for the single-particle 
spectrum of 2D liquid "^He has been found in the micro- 
scopic calculation of Ref. [2^ based on the method of cor- 
related basis functions^iiS^ (CBF).) According to Eq. H]), 
a spectrum e{p) having this shape implies that one can 
lower the energy of the system by introducing a quasihole 
at a certain momentum p < pp for which e{p) < 0, and 
a quasiparticle at a certain momentum p > pp for which 
eip) > 0. 

In the general framework of Landau theory, the quasi- 
particle spectrum is a functional 



SE 
5n{p) 



^J■ 



(3) 



of the quasiparticle momentum distribution n{p); hence 
the two quantities should be evaluated self-consistently. 
The functional approach as formulated and applied in 
Refs. 18,19,20 cannot be used in the 2D electron gas be- 
yond the critical point , since this approach has so far 
only been developed for states of Fermi systems describ- 
able in terms of standard Landau quasiparticles having 
the Fermi-step momentum distribution n-pi^{p). More- 
over, extension of this approach to finite temperatures 
remains incomplete. 

Nevertheless, an effective tool is available for self- 
consistent evaluation of the single-particle spectrum e{p) 
and momentum distribution n{p) beyond the point of in- 
stability of conventional Landau state, and indeed within 
Landau theory itself. This tool is the basic Landau for- 
mula 
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dp 
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dn{p') 
dp' 



dv' 



(4) 



relating the group velocity of quasiparticles, de/dp, to 
their momentum distribution 



n{p) 



1 



<P)/T 



1) 



(5) 



in terms of the Landau interaction function /(p, p')- 
The quasiparticle interaction function is the second vari- 
ational derivative of the energy functional E[n], 



/(p,p',W) 



6^E 



Sn{p) 5n{p') ' 



(6) 



and hence itself depends on the quasiparticle momentum 
distribution n(p). If non-zero temperatures are consid- 
ered, the T-dependence of this function should in prin- 
ciple also be taken into account. The theory presented 
here rests on two main assumptions. The first is that the 
difference between the interaction /(p, p',rs, [n]) corre- 
sponding to the true quasiparticle distribution n{p) at 
Ts > r^, and its counterpart /(p,p',rs, [tifl]) evaluated 
with Landau distribution nphip), is small as long as the 
relative size of the region of momentum space in which 
riFhip) is altered, remains small. Thus, we restrict our 
analysis to values of adjacent to the QCP. The sec- 
ond assumption is that the T-dependent terms in the 
quasiparticle interaction function produce contributions 
to the single-particle spectrum of order of T/e'p, where 
e'p = p\/2M . Thus, our analysis is further restricted to 
the regime of very low temperatures T/e^p <C 1. These 
two assumptions allow us to use the quasiparticle inter- 
action function obtained from the functional approach in 
the Landau relation (|4]), which becomes an equation for 
self-consistent evaluation of the quasiparticle spectrum 
and the quasiparticle momentum distribution. 

In Sec. II, we outline the approximation method em- 
ployed for evaluation of the quasiparticle interaction 
function of the 2D electron gas and the self-consistent 
scheme for subsequent calculation of e{jp) and n{p). 
Working within this framework, we then explore scenar- 
ios for rearrangement of single-particle degrees of freedom 
of the system just beyond its QCP (Sec. III). Section IV 
is devoted to calculation and analysis of the correspond- 
ing thermodynamic properties of the 2D electron gas near 
the QCP. In Section V, the quasiparticle phase diagram 
of the system is constructed and discussed. An Appendix 
presents further details on the calculation of the quasi- 
particle interaction function. 



II. SINGLE-PARTICLE SPECTRUM AND 
QUASIPARTICLE INTERACTION FUNCTION 

The well-known Feynman-Hellman formula, written as 



E^r-^ fde' fdq [Im x(g, u;) + 7rp<5(^)] (7) 











for the 2D electron gas, relates the energy E of the system 
to its response function xiQj <^)- This relation enables one 



to evaluate both the single-particle spectrum at T = 0, 
by calculating the first derivative of Eq. ([7]) as in Ref. [13, 
and the quasiparticle interaction function, as described 
herein. The latter is given by 



(52x(q,w) 













TT 



5n{p) Sn{p') 

u 

. (8) 

in term of the second derivative of x(cii ^) with respect 
to the occupation numbers. 

If the Landau quasiparticle interaction function is 
available, Eq. (|4]) can be treated as an equation for n{p). 
Specifically, making use of Eq. (O the Landau relation 
can be recast in the form 



T 



dn{p) p 
i{p)[l-n{p)] dp " M 



, dn{p') I 
fiP,P)-g^dv 



(9) 

with n{p) = — Jp {dn{p' ) I dp') dp' . 

A microscopic prescription for evaluation of the quasi- 
particle interaction function is given in the Appendix. Al- 
though full numerical realization of this scheme presents 
a formidable challenge, calculation of harmonics of the 
interaction function is tractable. Based on analysis of 
different contributions to /(p, p', r^, [jifl]) and computa- 
tion of its zeroth and first harmonics at the Fermi surface 
as functions of r^, we have developed a method for ap- 
proximate determination of the quasiparticle interaction 
function. To elucidate the essence of this approximation, 
we note first that if the momenta p and p' lie exactly 
at the Fermi surface, the interaction function depends 
only on the difference q = |p — p'|. Analysis has shown 
that the main contribution to this function has a strongly 
pronounced peak at qo ~ 1.97pF. Since our considera- 
tions are restricted to the immediate vicinity of the QCP, 
only momenta p and p' close to pp enter Eq. ([9]). Ac- 
cordingly, we approximate F(p, p') = iVo/(p,p'), where 
A^o = M/t:, in terms of the analytical form 



F{q) 



(10) 



which by construction is peaked at qo. To choose the 
other parameters specifying this function, we have made 
ab initio calculations of the zeroth and first harmonics 
of the interaction function and their derivatives with re- 
spect to Ts, and used the results, namely Fo{r^) = —1.3, 
Fi(rf^) = 1.0, and {rf / Fi{rf))dFi/ dr, = 1.05, as data 
to be fit by adjusting 7 and /3 in Eq. ([TU)) . In this 
way these parameters were determined to be 7(rs) = 
0.08(1 -I- imrs/rf) and [3 = 0.135. 

With the interaction function dependent only on |p — 
p'l, Eq. ([U reduces to the relation 

<P) = eo(p) + / /(|P - p'l) n{p') dv' , (11) 



in which £o{p) = p^/2M — /i is the bare sp spectrum. 
Eq. pTjl . together with Eq. ^ and the normalization 




FIG. 2: Single-particle spectrum e(p) in units of 10^^ e% (top 
panels) and derivative d€{p)/dp in units of v% = pp/M (bot- 
tom panels) for the 2D electron gas at rs = 6.8 (left column), 
rs = 6.9 (middle column), and = 7.0 (right column). The 
spectrum and its derivative are shown as functions oi p/pp at 
four (line-coded) temperatures expressed in units of e%-. 



condition 



n{p) dv — p , 



(12) 



provide a set of equations for self-consistent evaluation 
of the quasiparticle spectrum e(p) and momentum dis- 
tribution n[p). The chemical potential which enters 
Eq. pT|) as a constant of integration, is also determined 
self-consistently. 



III. REARRANGEMENT OF 
SINGLE-PARTICLE DEGREES OF FREEDOM 

We have solved Eq. (fTT|) for rs values situated in the 
vicinity of rf and on both sides of the QCP. To ana- 
lyze the nature of the solutions, we start with the results 
shown in Fig. 2 at r^. < r^, i.e., for densities above that 
of the QCP. As seen in the bottom panels of this figure, 
where the derivative de/dp is plotted against momentum 
p, the effective mass grows rapidly as approaches the 
QCP. The ratio M* /M at low temperature is about 6 
at Ts = 6.8, reaches 15 at = 6.9, and diverges at 
Ts = 7.0. We observe that near the QCP, the effective 
mass M*{rs) decreases with increasing temperature, and 
the closer is to r^, the more pronounced is the de- 
crease of M* (T) . The results shown in the bottom- right 
panel of Fig. 2 are in accord with the behavior of the 
group velocity predicted bythe inflection-point model of 
the QCP proposed in Rei.^, in that pf/M*[T) is found 
to vary approximately as T^/"^. 

Results from evaluation of the quasiparticle momen- 
tum distribution n{p) and spectrum e(p) at > rf, 
i.e., for densities below that of the QCP, are displayed in 
Figs. 3 and 4. Inspection of Fig. 3 informs us that at low 
temperatures in this region beyond the critical point, the 
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FIG. 3: Occupation number n{p) (top panels), single-particle 
spectrum e(p) in units of e% (middle panels), and ratio e(j>)/T 
(bottom panels) for the 2D electron gas at Vs = 7.1 (left 
column) and = 7.2 (right column). All three quantities 
are shown as functions oi p/pp at different temperatures ex- 
pressed in units of e'^. 

quasiparticle degrees of freedom are rearranged in a sce- 
nario involving the formation of a multi-connected Fermi 
surface. To the authors' knowledge, Frohlicb^ was the 
first to call attention to the possibility of this type of 
Fermi surface (see also Ref . Is^) . At Ts = 7.1 and = 7.2, 
the rearranged quasiparticle momentum distribution of 
the 2D electron gas at T = 0, for either spin projection, is 
defined by a fully occupied circle {n{p) — 1), surrounded 
by an empty inner ring {n{p) ~ 0), surrounded in turn by 
a fully occupied outer ring, all centered on p = 0. At the 
small finite temperatures indicated in Fig. 3, increasing 
T leads to population of the inner ring {n{p) > 0) at the 
expense of depletion of the outer ring (n(p) < 1). The 
spectrum appears to have a cubic shape 

with three zeroes. In this expression, the coefficient 
D{rs), being negative beyond the QCP, vanishes at = 
r^, while the coefficient A of the cubic term does not 
depend on . As the system moves away from the QCP 
at Vs > r^, the relative size of the momentum region 
affected by the quasiparticle rearrangement, measured 
V — {Pf~Pi)/PF, increases as 77 cx (r^/r^ — 1)-'^/^. 
With rising temperature, the momentum distribution 
Ti{p,T) and the ratio e{p,T)/T change dramatically. At 
T/e°p - 3 X lO-'^ (for = 7.1) and T/e^, - 10"^ (for 
Ts = 7.2) the values of n{p) in inner and outer rings of 
the T = distribution equalize, and the extrema in the 
spectrum e{p) disappear. 

As the temperature continues to increase, the quasi- 



particle spectrum and the momentum distribution each 
undergo a metamorphosis illustrated in Fig. 4. Inspec- 
tion of the temperature dependence of these quantities 
reveals a crossover from a state having a multi-connected 
momentum distribution created by formation of a ring- 
shaped "bubble," to a state characterized by a region of 
momentum space in which the ratio e{p,T)/T and the 
distribution n{p) are independent of temperature. The 
existence of such a finite domain over which temperature 
independence prevails is a signature of the phenomenon 
of fermion condensation, i.e., the formation of a fermion 
condensate (FC) in the momentum region involvedJ^i^i 
Comparing Figs. 3 and 4, it is seen that the size and 
location of the region of momentum space occupied by 
the FC coincide respectively with the size and location 
of the region in which the "bubble" rearrangement sce- 
nario - depopulation of the inner ring and population of 
the outer ring - operates at lower temperatures. 

The emergence, at low but not the lowest tempera- 
tures, of a sizeable domain in momentum space within 
which the occupation numbers differ from 1 and 0, im- 
plies that there is a crossover from a state with small en- 
tropy to another state with sufficiently large entropy. To 
some extent the phenomenon that occurs is analogous the 
Kosterliz-Thouless transition. Although we are deal- 
ing with a crossover rather than a real phase transition, 
this analogy enables us to estimate the temperature of 
the crossover. To do so, we compare the free energies 
per particle F = E ~ TS in the bubble phase and the 
FC phase. The entropy associated with the bubble rear- 
rangement is taken to be zero, while that corresponding 
to the FC is estimated as 77 = {pf—pi)/pF- The loss in 
energy due to the transition from the bubble phase to the 
FC phase is estimated^ as AE ~ \D{rs)\ri'^ Sp, and the 
coefficient of the linear term in Eq. (fT3)l . as |I?(rs)| ~ 
Thus, AE ~ Ty^e^. Since the difference AS* is of order 77, 
the free energy of the FC state becomes lower than that 
of the bubble state at Tq - AE / AS ~ \D{rs)\-qe%. ~ 
\D{r,)\^/^e%. 

Proceeding to higher temperatures, the region belong- 
ing to the FC gradually shrinks and finally disappears. 
To analyze the type of crossover taking place in this case, 
it is useful to study thermodynamic characteristics of the 
2D electron gas. 



IV. LOW-TEMPERATURE THERMODYNAMIC 
PROPERTIES OF 2D ELECTRON GAS 

The density of states 

NiT) = -j'-^d. (14) 

is a property well suited for analysis of the behavior of 
the 2D electron gas in different temperature regimes. The 
ratio of the inverse of this quantity to its inverse for the 
corresponding 2D Fermi gas, -/Vq"^ = tt/M, is plotted 
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FIG. 4: Same as in Fig. 3 but at higher temperatures. 



on a log-log scale in Fig. 5. The horizontal line corre- 
sponds to FL behavior, while the straight line with slope 
(or critical index) a corresponds to the NFL regime with 
N{T) cx T"". As seen in the top panel of Fig. 5, the 
FL behavior of the density of states at < persists 
only up to some characteristic temperature T*. In the 
region around T*, a crossover occurs from FL to NFL 
behavior, with a critical index a ~ 0.6. This result is 
in essential agreement with the inflection-point model of 
the QCP developed in Ref. [HI, which yields a = 2/3. 
The difference between 0.6 and 2/3 could be due to the 
term in the sp spectrum quadratic in p — pp, which is 
absent in the inflection-point model^^i According to the 
top panel of Fig. 5, T* - IQ-^eO, at r, = 6.8. As 
approaches the point = 7.0, the temperature T* de- 
creases and vanishes at the QCP. This is apparent in the 
middle panel of Fig. 5: at rg = 7.0 the inverse of the 
density of states, N^-^{T), shows universal NFL behav- 
ior with a critical index a ~ 0.6 within the entire range 
of temperatures under consideration. As demonstrated 
in the bottom panel of this figure, at > the quan- 
tity N~^(T) undergoes two crossovers separating three 
different temperature regimes. At low temperatures, it 
exhibits FL behavior with a large effective mass. In the 



vicinity of Tq (which is roughly 10 Sp at 



7.2), 



a crossover occurs from the FL regime to an intermedi- 
ate NFL regime. With further increase of temperature, 
another crossover occurs in the vicinity of T* (which is 
roughly 5 x 10"^ £^ at = 7.2), leading to the NFL 
regime with critical index a ~ 0.6. 

Returning to Figs. 3 and 4, we may clarify the rela- 
tionship of the three temperature regimes exhibited by 
the density of states, to the different regimes of behavior 
found in the quasiparticle momentum distribution n(j)) 
and spectrum e{p). The bubble phase, with its cubic-like 
spectrum (jl3p . evidently corresponds to the FL behav- 
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FIG. 5: Inverse of density of states N~^{T) in units of A'^,^^ 
for the 2D electron gas at Vs — 6.8 (top panel), Vs — 7.0 
(middle panel), and — 7.2 (bottom panel). Dashed lines 
show different temperature regimes. 



ior of N^^{T). Smoothing of the bubble and flattening 
of the quasiparticle spectrum e{p) in the region 
corresponds to the crossover to the FC state. However, 
the spectrum e(p) is also flat in the regions adjacent to 
the points pi and p f , and hence these regions also play a 
significant role in producing NFL behavior of thermody- 
namic properties. Indeed, the density of states. 



N{T) 



No[yf{T) + Un{T)] 



T 



decomposes into a sum of the FC term 

Pf 

Vf{T) = — / n^p) [1 - n*(p)] dp . 
PF J 

Pi 

and non-condensate part 

MT)^v>{T) + v<{T) 
consisting of two terms. Upon introducing 



i^„(r;ei,e2) 



f 

PF 



i{e) [1 - n(e)] de 



(15) 



(16) 



(17) 



(18) 



one can rewrite these two terms as — i^„(r; — /i, e^) 
and = i/„(T; e/, oo). At 77 ^ f the group velocity 
near the boundary pf can be approximated using 



dt{p pf)/dp ~ Vf{T) -f v^'^\p ~ PfY 



(19) 



To estimate the value of Vf{T), i.e. the slope of the 
FC "plateau" in the spectrum e{p), we substitute the 
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FIG. 6: Entropy per electron for the 2D electron gas as a 
function of temperature in units of £% at three values of rs. 



FIG. 7: Ratio 3C/Tir'^No for 2D electron gas at three values 
of the parameter rs is shown against temperature T in units 
of e%. 



momentum-independent result n^:{p) for the momentum 
distribution in the region into the Fermi-Dirac 

formula ([5]) and obtain 



e{p,T>To)^Tln 



1 - 71^ (p) 

n^{p) 



Pi <P <Pf 



(20) 



As seen in the bottom panels of Fig. 4, the width e{pf) — 
t{pi) = ef — ei oi the FC zone at T > To appears to be 
of order T and almost independent of rj for r/,„in ~ 10~^. 
Thus, at ry > 7/min the FC group velocity is estimated as 



\ dp Jj, i]pF 



Pi <P <Pf , 



(21) 



which in turn yields Vf{T) ^ T. The formula for the 
group velocity outside the FC region but near the point pi 
is analogous to (fT9|) with Vi{T) ^ T. Using these results 
in algebraic manipulations similar to those performed in 
Ref. ul], we arrive at the relation 



N{T) = NoT-^[iyf + , 



ri/31 



const. 



(22) 



expressed in terms of the dimensionless temperature r — 

The FC is seen to play an insignificant role at small 77, 
say 77 < T^/'^, where the density of states behaves as 



N{T) (X 



(23) 



In the opposite case, 77 > t^^'^, the FC contribution dom- 
inates in Eq. (|22p . In our calculation, the intermediate 
situation 77 ^ t^I'^ occurs, so the behavior of the density 
of states in the region To < T < T* corresponds to a 
superposition of two terms in Eq. of the same order. 

The NFL difference Ax(T,p) = xiT,p) - Xfl(/o) be- 
tween the magnetic susceptibility and Pauli susceptibility 
is determined as 



.1/3 



1 + Go 



(24) 



where Go, the zeroth harmonic of the spin-spin quasi- 
particle interaction function, is independent of density 
in advance of the QCP.^ The following behaviors of the 
magnetic susceptibility derive from the above findings for 
the density of states. At r > r* oc ry^, one has 



(25) 



whereas at t < r* the magnetic susceptibility of the 2D 
electron gas. 



X(T) « T-i , 



(26) 



imitates that of a gas of localized spins. At r < tq < t* 
the Curie-like behavior ((26)) is replaced by independence 
of T. 

Fig. 6 shows the entropy per electron on a log-log scale 
at three values of r^. Three different regimes may be 
discerned. Two of them are seen below the QCP, specif- 
ically at Ts — 6.8, namely (i) the FL domain with the 
usual linear temperature dependence and (ii) the NFL 
regime where S{T) oc T^"" with a ~ 0.6. At the critical 
point, i.e., at — 7.0, the slope of the curve S{T) is 
almost independent of temperature. Beyond the QCP, 
i.e., for rs > 7.0, all three regimes are present, (i) the FL 
regime at low T, (ii) the intermediate NFL regime cor- 
responding the emergence of the FC, and (iii) the NFL 
regime with a ~ 0.6. Similar observations follow from 
analysis of Fig. 7, where the ratio C{T)/T of the specific 
heat to temperature is depicted. The sharp drop of this 
ratio at rg — 7.2 corresponds to the crossover to the FC 
state, where the temperature dependence of the entropy 
is weak. 

The evolution of the temperature regimes with passage 
through the critical point are clearly revealed when the 
reduced Sommerfield- Wilson ratio 



3plC{T) 



(27) 



where Xo{T) = (1 -I- Go)x{T), is plotted against temper- 
ature for different rs values in the critical region. This 
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FIG. 8: Temperature dependence of the Sommerfield- Wilson 
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FIG. 9: Density of states N{rs) in units of A'^o for the 2D 
electron gas at four values of temperature T in units of £%. 



quantity is unity for the FL, and is usually strongly en- 
hanced in case of NFL behavior.— At rs = 6.8, the 
crossover from FL to NFL behavior is reflected in an 
enhancement of i?sw(r) for T > lO^e^p. At = 7.2, 
the ratio (P?]) goes to unity in the low-temperature limit, 
confirming the FL character of the bubble phase, whereas 
its enhancement to a value 7 and subsequent de- 
crease to ~ 3 with increasing temperature reflect the two 
crossovers discussed above. 

Up to this point, we have changed temperature at fixed 
values of Vg . Let us now go across the critical value at 
fixed temperature. The results for the density of states 
versus Vg evaluated at different temperatures are plotted 
in Fig. 9. Two different types of a behavior of N{rs) 
can be distinguished. Consider first the lowest of the 
four curves, corresponding to t = T/sp = 10^^. This 
curve shows a monotonic increase of the density of states, 
with a change of slope near the QCP, behavior which 
can be understood as follows. Below the QCP, i.e. at 
Ts < r^, the density of states N{T, r^) is proportional to 
the effective mass M*(T, rg) which, in the vicinity of the 
QCP, behaves as^i 

M*{rs,T) ^ 1 

Here, the positive constant D{rs) cx 1 — rg/r^ is the 
coefficient of the linear term (in p — pp) of the formula 
(|13|) for e(p), which provides a good approximation to 
the quasiparticle spectrum. In Eq. pSj) . v = (9A/4)^/'^ 
is determined by the rg-independent parameter A con- 
trolling the cubic term in Eq. (see Ref. '2T'), while k 
is a constant of order one. Beyond the QCP, a fermion 
condensate is present at T/e^p = 10^^ . The behavior of 
N{rs,T) in the state with a FC results from the inter- 
play of two contributions, as described by Eq. (P^ . The 
first contribution, associated with the FC, behaves like 
r]T~^ (x \D{rs)\^^'^T~^ , while the second, coming from re- 



gions adjacent to the FC domain, is proportional to t^^/'^ 
and independent of rg. 

We turn now to the uppermost curve in Fig. 9, corre- 
sponding to T = 3 X 10~^. Below the QCP, the behavior 
of this curve can also be understood in terms of Eq. (|28p . 
Beyond the QCP, the bubble rearrangement is favored 
by the low temperature, and the resulting FL behavior of 
the density of states can be again described by Eq. (|28p , 
with all three sheets of the Fermi surface taking part and 
making the coefficient k considerably larger than unity. 

An essential message of Fig. 9 is that at sufficiently low 
temperature, the effective mass of the 2D electron gas ex- 
tracted from an analysis of its thermodynamic properties 
should decrease beyond the QCP. 



V. QUASIPARTICLE PHASE DIAGRAM 

The information we have gained on the behavior of 
quasiparticle degrees of freedom of the 2D electron gas 
near the quantum critical point may be integrated into a 
quasiparticle phase diagram. We do this with the caveat 
that our treatment has not taken into account the feed- 
back of the modification of single-particle degrees of free- 
dom on any pertinent collective modes. Accordingly, 
what we obtain is only a "quasiparticle sketch" of the 
real phase diagram of the system. 

Phase diagrams of other strongly correlated systems 
are customarily drawn in a T — x plane, where x is some 
external parameters, notably pressure, doping, or mag- 
netic field. Instead of these, we have considered varia- 
tions of the spacing parameter rg measuring the extent 
to which the system is rarefied. Thus, in Fig. 10 the 
phase diagram of the 2D electron gas is drawn in (T, x) 
variables, with x — 1 ~ rg/r^ . On the top axis of the 
plot, we show also the coefficient D{x) oc a; of the lin- 
ear term in Eq. (|13p for the quasiparticle spectrum. The 
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FIG. 10: Quasiparticle phase diagram of the 2D electron 
gas in (r, x) variables with x — 1 — rs/r'^. The dashed 
curves show crossovers between the usual Fermi liquid (FL), 
the Fermi liquid with bubble (FL-B), the fermion condensate 
(FC) state, and the non-Fermi-liquid (NFL) phase with crit- 
ical index a ~ 0.6. The parameter D[x) is indicated on the 
top axis. 



latter parameter is positive below the QCP and negative 
above; its vanishing at the QCP reflects the emergence of 
the horizontal inflection point in the quasiparticle energy. 
Beyond the QCP, the negative sign of D{rs) corresponds 
to three roots of the equation e{p) — 0, and the presence 
of a bubble in the quasiparticle momentum distribution. 

Inspecting the phase diagram presented in Fig. 10, we 
start from the side to the right of Xoo, i-e., below the 
QCP, where D{rs) > 0. Here the low-temperature re- 
gion of the diagram is occupied by an ordinary Landau 
Fermi liquid. The foregoing analysis tells us that at a 
temperature around T*, a the crossover occurs from the 
FL regime to a universal NFL domain that occupies the 
upper part of the diagram. The dependence T*(D) can 
be estimated using Eq. As a function of T, the 

behavior of this expression changes at a temperature for 
which vt'^/'^ ^ D. Thus we estimate that the tempera- 
ture of the crossover to the NFL behavior in the domain 
X > goes like 



T*{D>0)^:^-^/^\D\^/^e°p 



(29) 



This temperature is traced by the dashed line in the right 
side of the phase diagram of Fig. 10, separating the light 
gray FL region from the white NFL region. 

Now consider the left side of the phase diagram lying 
beyond the QCP, where D{x) < and two crossovers 
occur. Under increasing temperature there is first a 
crossover from the FL state having a bubble (FL-B) to 
the FC state. An estimate of the temperature Tq of this 
crossover was made in Sec. Ill, with the result tq ~ \D\ri. 
The dependence of 77 on can be easily found by solving 
the equation e{p) = 0, which identifies the boundaries of 
the region in which the Landau quasiparticle distribution 
is altered. For the spectrum as represented by Eq. ([13]), 



this equation has the explicit form 
Dt] + v'^rf/'d = 

yielding 

r/ = 3i.-3/2|^|i/2 



(30) 
(31) 



Incorporating this dependence into tq ^ \D\'q^ we obtain 
To{D < 0) - 3i/-3/2|^|3/2 (32) 

The crossover in question is indicated in Fig. 10 by the 
curve in the left part of the phase diagram separating the 
light-gray FL-B regime from the gray region correspond- 
ing to the FC state. 

The second crossover, from the FC regime to the uni- 
versal NFL regime with critical index 0.6, occurs as the 
temperature increases up to the value T*{D < 0) at 
which the FC contribution to thermodynamic character- 
istics of the system and that of regions adjacent to the 
FC territory become comparable. For example, the FC 
contribution to the density of states has been estimated 
in Sec. IV as r]T~^, and that from adjacent regions, as 



j,-i^-2/3_ rpj^jg implies that t*{D < 0) 
hence 



ly^ij'^ and 



T*{D < 0) - 27 J/" 



3/2 1^,3/2^0 



(33) 



The resulting dependence of T*(D < 0) on x is shown 
by the upper curve in the left part of the phase diagram, 
which separates the FC region from the NFL region. 

As mentioned at the beginning of this section, the 
phase diagram we have generated is a semi-finished prod- 
uct. The next step in building a complete phase diagram 
of the 2D electron gas is to characterize the feedback from 
rearrangement of the single-particle degrees of freedom 
on the collective modes. It is expected that one of the 
Pomeranchuk stability conditions is eventually violated, 
leading to collapse of the corresponding collective mode. 
Preliminary analysis indicates that this violation may be 
due to unusual behavior of the quasiparticle polarization 
operator beyond the QCP. Future work will aim at deter- 
mining which channel of instability is the first to appear. 



VI. CONCLUSIONS 

We have studied low-temperature properties of the 2D 
electron gas at densities close to the point p^o of diver- 
gence of the quasiparticle effective mass. A number of 
significant results have been obtained by means of a mi- 
croscopic approach based on fundamental relations of the 
Landau Fermi-liquid theory and an approximate evalua- 
tion of the quasiparticle interaction function. 

For Ts < rf, i.e., for values of the dilution parameter 
below the point of divergence of the effective mass, we 
find that the temperature T* up to which the 2D elec- 
tron gas continues to exhibit FL thermodynamic prop- 
erties, decreases with increasing and vanishes as 
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approaches r^. At T ~ T*, a crossover occurs to uni- 
versal NFL behavior of thermodynamic properties: the 
density of states N{T), the spin susceptibihty x(r), and 
the ratio C{T)/T of the specific heat to temperature scale 
universally as with critical index a ~ 0.6. This find- 
ing is consistent with predictions of the inflection-point 
model of the quantum critical point (QCP).^"'^ 

We have evaluated the spectrum e{p) and the momen- 
tum distribution n{p) of quasiparticles beyond the QCP, 
i.e. at Ts > rf^- In this regime, the single-particle de- 
grees of freedom are altered, more or less profoundly, 
relative to the standard FL picture. At very low tem- 
peratures, the new quasiparticle momentum distribution 
is multi-connected. An empty ring appears in the oc- 
cupied Fermi circle, as the self-consistent quasiparticle 
spectrum e{p) acquires a cubic-like shape with three ze- 
roes at three boundary momenta. This state, however, 
possesses the usual Fermi liquid thermodynamic prop- 
erties, in that N{T) is temperature independent, as are 
X(T) and the ratio C{T)/T. Moreover, the Sommerfield- 
Wilson ratio is unity, as it should be in a FL state. 

When the temperature rises to a certain value Tq, a 
crossover to a state having a fermion condensate is found 
to occur. This more radical rearrangement of the single- 
particle degrees of freedom is distinguished by fractional 
occupation of momentum states, hence essential depar- 
ture from the conventional Landau Fermi-liquid picture, 
over a finite range of momenta where the fermion conden- 
sate resides. The momentum distribution T>To) in 
this domain, as well as the ratio e{p,T)/T, proves to be 
independent of temperature; the spectrum €{p, T) itself, 
being proportional to T, proves be almost flat. The tem- 
perature dependence of the density of states N{T), and 
likewise for the spin susceptibility x(T'), is determined 
by an interplay of a pure FC term proportional to 
and a NFL term proportional to T^'^^'^. The tempera- 
ture dependence of the entropy is weak, while the specific 
heat decreases sharply and the Sommerfield- Wilson ratio 
jumps up to about 7. With further increase of tempera- 
ture, a second crossover is found, going from the FC state 
to the state with universal NFL behavior. 

Assembling results from the analysis of thermody- 
namic properties, we have sketched a quasiparticle phase 
diagram of the dilute 2D electron gas in the immediate 
vicinity of the QCP. There are four regions in this dia- 
gram: (i) a conventional one corresponding the FL be- 
havior of ordinary Landau quasiparticles, (ii) the region 
associated with the multi-connected quasiparticle distri- 
bution where FL temperature behavior persists, (iii) the 
domain created by the FC rearrangement, and (iv) the 
regime of universal NFL behavior. 

Also of interest is our finding that at a fixed temper- 
ature for which the multi-connected quasiparticle distri- 
bution is favored, there is an enormous enhancement of 
the effective mass M* as increases to the QCP, but 
no true divergence. Moreover, beyond r^, this enhance- 
ment is followed by a decrease of M* . At temperatures 



for which the FC state is favored, the rate of growth of 
the effective mass decreases with increasing beyond 
the QCP. 

New experiments using ultraclean Si-MOS samples 
could clarify the nature of the QCP in 2D electronic sys- 
tems inhabiting these structures. 
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APPENDIX 

In this appendix we describe the scheme adopted for 
evaluation of the quasiparticle interaction function as ex- 
pressed in Eq. ([5]) . A microscopic functional approach to 
this problem begins with the connection 



1 - R{q,uj)xo{q,uj) 



(34) 



of the linear response function x{1^ ^) entering the in- 
tegrand in the r.h.s. of Eq. ([8]) to the response func- 
tion Xol^j ^) of the noninteracting gas, through the effec- 
tive interaction R(q,u;). Accordingly, variation of the re- 
sponse function xiQ, ^) affects both the function Xo(<Z: ^) 
and the effective interaction R{q,uj): 



<5^x(a,w) 



5n{p) 5n{p') (5n(p) (5n(p') 



-2ip^{q,uj)R{q,uj[ 



Sn{p) Sn{p') 



terms containing 



SR{q,uj) 
Sn{p) 



(35) 



where ^p{q,uj) = [1 - i?(g, t^)xo(g, w]! ^ 

The analysis performed in Ref. |20 indicates that the 
shape of the effective interaction has a rather smooth 
dependence on r^. For this reason, the terms contain- 
ing the variational derivative SR{q,Lu)/Sn{p) serve only 
to produce a slight renormalization of the chemical po- 
tential, without contributing significantly to the shape of 
the spectrum. By the same token, these terms are also 
neglected in the present treatment. 

We consider first the contribution /^^•'(p, p') to the 
interaction function resulting from the product of first 
variational derivatives of the fimction xo('Z;'^): namely 
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(1 - npL(p-q)) (1 - npL(p'-q)) <5(eO - e° - e° ^ + _ J ^'{q, 6° -."J 



(l-nFL(p-q))(e°-e°_e 



(e°- 



(l-nFL(p'-q))(e^,-60,_, 



^ p 



p-q 



2^ K-e'r+w^] [(e 



■ <y5'^(q, iw) 



(36) 



Following the strategy discussed in Section II, we as- 
sume that both of the momenta p and p' are located 
near the Fermi surface. We also assume both of them 
to be larger than pp^ as is conventionally done in other 
problems of this kind. Thus, we deal with the function 
f^^\x) = f^^\pFTPF,x), where x = cos(p,p'), and we 
are interested in the region a; ~ — 1 that corresponds to 
backward scattering and provides the main contribution 
to the shape of the spectrum near the Fermi surface. It 
should be kept in mind that within the diluteness regime 
under consideration, i.e., ~ r^, the function Lp{q,{)) 
possesses a sharp maximum at q ~ 2pi?, such that the the 
dominant contribution to the integral in Eq. (j36p comes 
from the region where q ^ 2 pp. Assuming — 6^, — 0, 
we note that the first term of the integrand of Eq. 
contains the delta-function (5(ep_q — ep'_q)- For x ^ —1, 
the condition of equality of the energies fp-q and ep/_q 
suppresses the integral over angles of the vector q, since 
it receives contributions only from directions of q tan- 
gent to the Fermi surface at the points p and p'. The 
first term of the integrand in Eq. (|36)) yields appreciable 
contributions to /^^•'(x) only at x ~ 1, noting that in this 
case all directions of q carry similar weight. The second 
and third terms of the integrand contain expressions 



(e^q)'-(4'-q)' 



and 



p'-q 



y p 



p q 



respectively, thus implying that their contributions peak 



at ep_q — fp'-q- The same considerations as for the 
first term show that neither of these terms contributes 
noticeably to f^^\x) at x ~ —1. In the fourth term, 
the region w ^ 1 dominates the integration over the 
variable w. When evaluated within the approximation 
i?((7, lu^l) — R{q,Q) (which overestimates the integral), 
the fourth term produces a contribution 



^0 



^^.-q + ^^ 



(37) 



p^-q 



sponding contribution to /(.t~ — 1), made by replacement 
of the denominator in Eq. (|37p by its approximate value 
q^/A/, indicates that the effect is rather small, namely 
about 10 — 15% at ~ 7. This contribution has not 
been neglected, but instead is evaluated approximately. 

We turn now to evaluation of the main term in the 
interaction function, denoted /^^HPiPOj containing the 
second variational derivative of the response xo(9,^)- 
This derivative is easily evaluated as 

1 5'^X{<1,^) 



2 2(5n(p) 5n(p') 



and we then obtain 



e°,-c.)<5(q-p+p') , 

(38) 



f^'Hp,p')^-J de^^^^^(p-p',e°-e°,) . (39) 



For momenta p and p' lying at the Fermi surface, 
Eq. ([55]) reduces to the expression 



2n 

q 



(40) 



the integral over being evaluated numerically. 

The interaction function f''^\q) calculated for the 2D 
electron gas at = 6, 7, and 8 is displayed in Fig. 11. 
The absolute value of this function is seen to possesses a 
strongly pronounced maximum at q = gc — 2pi?, which 
increases with increasing r^. The peak and its enhance- 
ment with increasing are related to the proximity of 
the system to an instability against spontaneous appear- 
ance of ch arge-d ensit y waves (CDW) with wave vector qc 
(see Refs. l20li26ll27ll28l ). This situation is known to be re- 
sponsible for divergence of the effective mass.^^'^^ In this 
regard, it is worth noting that the negative sign obtained 
for the quasiparticle interaction function f(q^2pp) pro- 
vides the correct positive sign of the first harmonic. 



F,=No fix) 



to the interaction function. An estimate of the corre- 



E dx 



(41) 
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which equals 1 at the QCP. At the same time, the nega- 
tive sign of the quasiparticle amphtude does not support 
the idea that the quasiparticle weight Z vanishes near a 
CDW instability. The latter eventuality is possible only 
in case the quasiparticle amplitude, divergent at q ~ 2pf , 
has a positive sign, as erroneously assumed in Ref. |l6| in 
considering enhancement of the amplitude in the vicinity 
of such an instability. 



1 2 3 

FIG. 11: Quasiparticle interaction function, calculated within 
the microscopic functional approach for the 2D electron gas 
and plotted as f{q) No at three values of the diluteness pa- 
rameter Ts- 
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